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We report on polariton condensation in photonic molecules formed by two coupled micropillars. 
We show that the condensation process is strongly affected by the interaction with the cloud of 
uncondensed excitons. Depending on the spatial position of these excitons within the molecule, 
condensation can be triggered on both binding and anti-binding polariton states of the molecule, on 
a metastable state or a total transfer of the condensate into one of the micropillars can be obtained. 
Our results highlight the crucial role played by relaxation kinetics in the condensation process. 

PACS numbers: 71.36.+C, 67.85.Hj, 78.67.Pt, 78.55. Cr 



Most of the experimental studies in atomic Bose con- 
densates have explored conditions of thermodynamic 
equilibrium since typical condensate lifetimes are much 
longer than interaction times. Recent theoretical pro- 
posals have shown that out of equilibrium bosonic sys- 
tems present qualitatively new behaviors One pro- 
posed way to reach this regime is the use of photonic sys- 
tems with effective photon-photon interactions and dissi- 
pation provided by inherent optical losses @ ■ Localized 
to delocalized phase transitions 0, H| , highly entangled 
states |5[, or fermionisation effects in a ring of coupled 
sites [6J] are predicted in such systems. 

Microcavity polaritons are a model system for the in- 
vestigation of the physics of driven-dissipative boson con- 
densates They are the quasi-particules arising 
from the strong coupling between excitons confined in 
quantum wells and the optical mode of a microcavity. 
Because of their light-matter nature, polaritons present 
peculiar properties: they interact efficiently with their 
environment through their excitonic part [3, EH while 
their photonic part enables efficient coupling with the free 
space optical modes. Polariton condensates can be gen- 
erated in zero dimensional micropillars lll f or in arrays of 
pillars with fully controlled coupling [ly, ll7| ■ In this con- 
figuration, the non-equilibrium nature of polariton con- 
densates should allow the realization of metastable col- 
lective states, such as the self-trapped states in a bosonic 
Josephson junction fl8l - 420l | . 

In the present paper we investigate polariton conden- 
sation in photonic molecules obtained by coupling two 
micropillars. We demonstrate that polariton interactions 
strongly affect the way condensation occurs in such cou- 
pled system, not only modifying the wavefunction of the 
polariton condensate, but also the relaxation dynamics. 
This effect, specific to an out-of-equilibrium bosonic sys- 
tem, is illustrated by considering different positions of the 
non resonant excitation within the molecule. When the 
excitation spot is placed at the center of the molecule, 
polariton condensation is observed on both binding and 



anti-binding states. Interactions induce strong changes 
in the condensate wavefunction, the most important one 
being the change in its spatial anisotropy. 

When the excitation spot is positioned on one of the 
two coupled micropillars, condensation occurs in a very 
different way. As the excitation power is increased, the 
polariton condensate is first created in a metastable state 
localized in the excited pillar by an effect similar to self 
trapping fig. Further increasing the pumping power, 
the condensation tends to progressively evolve from a 
kinetic regime to a regime closer to thermodynamical 
equilibrium: massive occupation is observed on the low- 
est energy state of the system, mainly localized in the 
non-excited pillar. These results are simulated using 
a relaxation model including semi-classical Boltzmann 
equations and nonlinear Schrodingcr equation in a self- 
consistent way. The key features of our experiments can 
only be reproduced if one properly includes the changes 
of the relaxation rates due to interaction-induced modi- 
fications of the overlaps between the different polariton 
states. 



The microcavity sample, described in Rcf . ll| , consists 
in a A/2 cavity with a quality factor exceeding 16 000 and 
containing 12 GaAs quantum wells. Coupled micropillars 
were fabricated using electron beam lithography and In- 
ductively Coupled Plasma dry etching. The diameter d 
of the micropillars ranges from 3 to 4 fim and their cen- 
ter to center distance dec is varied from 2.3 to 3.7 /Ltm, 
corresponding to an expected coupling constant between 
0.1 and 1 meV. Note that dec is always kept smaller 
than twice the radius of the pillars, ensuring the direct 
coupling of the polariton modes of the two micropillars. 
A scanning electron micrography of an array of such pho- 
tonic molecules is shown in FigHJa. Microphotolumines- 
cence experiments are performed on single molecules us- 
ing a single mode cw TbSaph laser focused onto a 2 /im 
diameter spot with a microscope objective. The sample 
is maintained at 10 K and the excitation laser energy is 
tuned typically 100 meV above the lower polariton reso- 
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FIG. 1: (Color online) (a) Emission spectrum measured on 
a single molecule at low excitation power (d = 4fj,m and 
dec = 3.73 /-tin); A and B indicate the binding and anti- 
binding states. Inset : scanning electron micrograph of an 
array of pillars and molecules, b) Emission spectra measured 
on a 4jj,m round micropillar (black line) and on photonic 
molecules with d = 4^im and various values of dec- c) left 
(resp. right) column : measured (resp. calculated) emission 
pattern of the four lowest energy modes of a photonic molecule 
with d = 4 fim and dec = 3.73 \im\ d) Schematic of the 
hybridization of the individual pillar modes within a photonic 
molecule. 



nance, thus providing non-resonant optical excitation of 
the system. The emission is collected through the same 
objective and imaged on the entrance slit of a monochro- 
mator. The spectrally dispersed emission is detected with 
a nitrogen-cooled CCD camera. Wc define the detuning 
S = E c — E x as the energy difference between the lowest 
energy photonic mode and the exciton resonance. 

The polariton modes in these photonic molecules are 
investigated by photolumincscence measurements at low 
excitation power. An example of an emission spec- 
trum measured on a single molecule is presented in 
Fig. [TJa. Discrete emission peaks are observed corre- 
sponding to polariton quantum states fully confined in 
the microstructure. The broader line at higher energy 
is due to emission of the excitonic reservoir. The two 
lowest energy modes (labeled B and A) are attributed to 
the binding and anti-binding states arising from the hy- 
bridization of the lowest energy mode of each micropillar. 
The splitting between these states is proportional to the 
coupling between the two micropillars. It can be contin- 




FIG. 2: (Color online) Excitation at the center of 

the molecule: a)-c) Spectrally resolved emission distribution 
along the molecule axis measured for several excitation pow- 
ers; d)-f) Calculated emission patterns for several excitation 
powers, d = 4 \im and dec = 3.46 fim. S — —3meV 



uously tuned by changing the center to center distance 
dec, a s illustrated in Fig.[TJb. The left column of Fig.[TJc 
shows the spatial distribution (near field) of the four low- 
est energy modes in a molecule made of two 4 /xm pillars 
with dec = 3.73 fim (the first two panels correspond to 
the B and A modes). The right column shows the calcu- 
lated polariton wavefunctions considering a confinement 
potential taken as infinite outside the photonic molecule 
and equal to zero inside. As schematically illustrated in 
Fig. [T]d, all these states result from the hybridization of 
the optical modes of each individual micropillar, showing 
the strong analogy between our system and the orbitals 
of a diatomic molecule. 

We now discuss polariton condensation in these 
molecules under non resonant optical excitation. As we 
have shown in our recent work on single micropillars plj , 
such an excitation scheme not only populates the con- 
fined polariton states, but also creates a population of 
uncondensed excitons called the excitonic reservoir. Be- 
cause of the limited diffusion length of excitons, this 
reservoir remains localized in the excitation area. Re- 
pulsive interactions between polaritons and the excitonic 
reservoir strongly influence the precise quantum state in 
which polaritons accumulate. In the following we will 
consider two different locations of the excitonic reservoir 
within the molecule, which can be selected via the posi- 
tion of the excitation laser spot. First we will consider 
excitation conditions where the reservoir of uncondensed 
excitons is at the center of the molecule. Then we will 
address the case of asymmetric excitation, in which the 
excitonic reservoir is injected only in one of the two mi- 
cropillars. Specific spatial behavior of the polariton con- 
densates is observed in each case, driven by the polariton 
interaction with the excitonic reservoir. 

Figure [5]a-c present the measured emission distribu- 
tion along the axis of a photonic molecule excited at its 
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FIG. 3: (Color online) Excitation at the center of the 

molecule: a) Integrated intensity measured as a function of 
the excitation power; b) Splitting between the two lowest en- 
ergy emission lines measured as a function of the excitation 
power; c) Lowest energy emission lines measured with a polar- 
ization parallel (solid line) and perpendicular (squares) to the 
molecule axis for P = 0.05 P t h\ d) same as c) for P — 7 Pth- 
The inset shows the calculated spatial shape of the binding 
state 



center. As summarized in Fig. [3]a, above a well defined 
excitation threshold, a strong non-linear increase of the 
emission intensity is observed, one of the signatures of 
polariton condensation 
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Under this excitation con- 
dition, a massive accumulation of polaritons occurs both 
in the binding and anti-binding lowest energy polariton 
states. Interestingly, as the excitation power is increased, 
a progressive change of the spatial shape of the polari- 
ton wavefunction is observed, with vanishing probability 
density at the center of the molecule. This behavior is 
due to repulsive interactions with the excitonic reservoir 
injected in that region. 

To describe our experiments, we simulate polariton re- 
laxation using coupled semi-classical Boltzman equations 
and nonlinear Schrodinger equations in a self-consistent 
As described in the supplementary informa- 
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way 

tion, on one hand, a nonlinear Schrodinger equation is 
solved at each time iteration, in the presence of the po- 
tential induced by the excitonic reservoir and by the oc- 
cupation of the polariton states, to find the eigenstates 
of the system. On the other hand, the time evolution of 
these occupation numbers are calculated using the semi- 
classical Boltzmann equations in which the scattering 
rates are obtained using the exact shape of the polariton 
states obtained from the non-linear Schrodinger equa- 
tion. Both exciton-phonon and exciton-exciton scatter- 
ing mechanisms are considered. The reservoir is assumed 
to be thermalized at the lattice temperature. The key in- 
gredient to describe our experiments is to consider scat- 



tering rates proportional to the spatial overlap between 
initial and final states. For instance, the scattering rate 
of two excitons from the reservoir resulting in one po- 
lariton in state i and one exciton in the reservoir is pro- 
portional to J \tpi(x,y)\ 2 (p R (x,y)) 3 dxdy, where ipi(x,y) 
is the wavefunction of the polariton state i and pn{x,y) 
is the spatial distribution of the excitonic reservoir, de- 
scribed by a gaussian of width wr. To fit the data we 
use an exciton-exciton scattering rate Wxx = 2 10 3 s _1 
and an exciton-phonon scattering rate Wxp = 10 9 s _1 . 
The lifetime of the reservoir is equal to 400 ps, and that 
of the polariton states to 30 ps, a value extracted from 
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coherence measurements in ref. 

The results of the model for central pumping are shown 
in Fig. [2jl-f. The potential of the reservoir increases with 
the optical pumping and has a maximum at the center 
of the molecule, inducing the spatial separation of the 
states of each micropillar. As a result, both the cou- 
pling and the associated splitting between the two lowest 
energy states decrease. Indeed, this feature is shown in 
Fig. Ob: simultaneous to the spatial separation of the po- 
lariton distributions, the measured splitting goes down to 
a minimum at about P = 4 P t h ■ At higher powers, the 
measured splitting tends to increase again. In the follow- 
ing, we will show that at such high density the observed 
splitting is no longer related to the binding anti-binding 
splitting but is an anisotropic splitting induced by the 
lateral spatial shrinking of the condensate wavefunction 
in each micropillar. This feature is not reproduced by our 
model which does not include the polarization degree of 
freedom. 

Polarization resolved measurements were performed 
for different excitation densities. Figure [3]c and 03 d dis- 
play emission spectra linearly polarized parallel and per- 
pendicular to the molecule axis, below and above thresh- 
old. In the low density regime (Fig. |3lc), the presence 
of uncondensed excitons has a negligible effect, and we 
observe a polarization splitting in both the binding and 
anti-binding states. As previously reported theoretically 
and experimentally (23[, this splitting is larger for the 
binding state and amounts to 70/ieU. The inset of 
FiglHlc shows that the overall shape of the polariton bind- 
ing state is elongated along the molecule. This is why 
the lowest energy line of the binding state doublet is po- 
larized parallel to the molecule axis. The situation is 
different at high density (P > AP t h), where a stronger 
potential barrier is induced by the excitonic reservoir at 
the center of the molecule. In this case, the coupling be- 
tween the states of each micropillar becomes too small to 
be resolved in our experiment. The remaining splitting is 
then purely a polarization splitting. In each micropillar, 
the polariton wavefunction is shrunk by the excitonic re- 
pulsive potential and thus it is more elongated perpendic- 
ularly to the molecule axis. This change in anisotropy is 
evidenced in the polarization of the emission. The lowest 
energy emission line is now strongly polarized perpendic- 
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FIG. 4: (Color online) Asymmetric Excitation: a)- 

c) Spectrally resolved emission distribution along the molecule 
axis measured for several excitation powers; d)-f) Calculated 
emission patterns for several excitation powers, d = 3.5 fim 
and doc = 3.1/xra. 5 — —3meV 



ular to the molecule axis. 

Different dynamics is observed when the laser spot is 
positioned on top of one of the micropillars forming the 
molecule. Emission distributions measured under these 
conditions are reported in Fig. H]a-c. Close to thresh- 
old, condensation occurs in the pillar that is optically 
pumped. Interparticle interactions blue shift the energy 
of the ground state in the pumped pillar, decoupling it 
from the states of the other pillar , and severely limiting 
the Josephson transfer of particles from one pillar to the 
other. This localisation is very analogous to the origi- 
nal self-trapping effect, except that here the blue shift 
is mainly induced by the interaction between the exci- 
ton reservoir and the condensate and not by interactions 
within the condensate itself. As the excitation power 
is increased, the better relaxation kinetic and the pres- 
ence of excited states in the un-pumped pillar destroy 
the metastable state [24| and allow the system reaching 
its ground state. Condensation then occurs in the non- 
pumped micropillar (Fig. |4jc). 

We can reproduce this overall behavior using the self- 
consistent polariton relaxation model presented above. 
As shown in Fig. 2]d, & t l° w power the reservoir poten- 
tial creates an asymmetry in the two bottom states of the 
molecule, whereas the upper states arc not affected. The 
populations of the excitonic states in the reservoir are 
small and their relative distribution is given by the lattice 
temperature. When pumping increases (Fig. Hie), con- 
densation starts at the eigenstate possessing the largest 
overlap integral with the reservoir. This corresponds to 
the third polariton state (population n^). The reser- 
voir potential also significantly perturbs the two bottom 
states, which become both confined in the right pillar. 
This blocks the direct relaxation from the reservoir to- 



wards these states and further enforce condensation in 
the third metastable state localized in the pumped pil- 
lar. At high pumping (Fig. HJf ) a new effect is observed: 
condensation is triggered onto the ground state (of popu- 
lation ni). This does not arise from direct scattering from 
the reservoir, since the overlap between the ground state 
and the reservoir remains very small. This condensation 
is the result of efficient scattering from the intermediate 
third state toward the ground state via acoustic phonons. 
This regime occurs when the phonon assisted scattering 
rate Wxpn-i(\ +ni) from the intermediate state into the 
ground state overcomes the ground state radiative losses 
n\T\. This happens when 713 100 and further increase 
of the pumping does not change 713, but only increases ri\. 
The second quantized state remains weakly populated: 
it has a smaller overlap with the reservoir than the third 
state and a higher energy than the ground state. Over- 
all the complete condensation dynamics can only be un- 
derstood considering the modifications of the relaxation 
rates induced by interactions (as shown in the supple- 
mentary [22| when constant scattering rates are consid- 
ered, condensation in the metastable states can not be 
reproduced). 

To conclude, we have demonstrated polariton conden- 
sation in photonic molecules, a fully controlled system 
in which the coupling constant between two condensates 
can be adjusted. Depending on the precise location of 
the excitonic reservoir, strong renormalization of the po- 
lariton states is induced, resulting in strong modification 
of the condensation dynamics. In particular condensa- 
tion can occur in a metastable localized state , or in 
the ground state of the system. These results open the 
way towards a detailed investigation of Josephson oscil- 
lations (2o| or, more generally, the pumped-dissipative 
physics of arrays of coupled condensates in an engineered 
environment. 
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Abstract 

We describe in this supplementary information the model we have developed to simulate polariton 
relaxation in photonic molecules under the excitation conditions described in the paper. The results 
of these calculations are shown on figure 2 and 4 of the paper. Finally we show that the overall 
condensation dynamics is not reproduced if one does not take into account the interaction induced 
changes in the relaxation rates. 

PACS numbers: 
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We simulate polariton relaxation in a self-consistent way, using coupled semi-classical 
Boltzmann equations and nonlinear Schrodinger equations. At t — the system is empty 
and the pumping is switched on. We calculate the time evolution of the system until the 
establishment of a steady state. The nonlinear Schrodinger equation written below is used 
to find the spatial shapes and energies E itt . of the polariton wave functions ipi ttj {x,y). 

h 2 

+ap R (x,y)N R (tj) + a 

m is the polariton effective mass. This equation is solved at time step tj to find 
the eigenstates in presence of the potential created by the exciton reservoir (written as 
apn(x,y)Nji(tj)) and by the occupation of each polariton quantized states (last term of the 
equation), a = QxE^clb) 2 is the exciton-exciton interaction constant, where x the exciton 
fraction of the considered quantum state (1 in the reservoir, 0.5 in the polariton states), 
the exciton binding energy, ag the exciton Bohr radius. Pr(x, y) is the normalized reservoir 
density (J pn(x,y) dxdy = 1), N R (tj) the number of excitons in the reservoir. 

The reservoir is assumed to be thermalized at the lattice temperature: nn(E,tj) = 
n Fi{Q,tj) exp (—EftsT). The total number of particles in the reservoir is N R (tj) = 
2 ur(E, tj), where 2 stands for the spin degeneracy of bright heavy-hole excitons. Since 
excitons in the reservoir are confined in a quantum well, N R (tj) = nR(0,tj)SmkBT/2irh 2 , 
where S is the surface covered by the reservoir. Pr(x, y) is assumed to be a Gaussian function 
with a width wr close to the one of the excitation laser spot. 

The time evolution of rij, the occupation numbers of the quantized polariton states and of 
Nr, the reservoir population are calculated using semi-classical Boltzmann equations, taking 
into account both exciton-exciton and exciton-phonon interactions, as well as the presence 
of the thermalized excitonic reservoir. 

For the quantized polariton states i the equation reads: 



.tj/h 



3-1 



2 



lT = -T 1 - n i E (tfc WA^ + Wf P ) (n fc + 1 



+(rii + 1) E (WrxWWfl + W^pJ^n* 



1, fc < z 

e -E lk /k B T^ ^ > i 
e -E lk /k B T^ ^ < i 



l,k > i 

+ (WGcxiiflRH^ + Wf P I iR ) (N R (m + 1) - ni^Ie-WfcBf) 
and for the reservoir state R the equation is: 



dN 



R 



dt 



(w xx i iRRR N R + W f P i iR ) (n r ( ni + 1) - y -wj 



(2) 



P is the non-resonant pumping populating the reservoir, rp = 400 (resp. Tj = 30 ps 
) is the excitonic reservoir (resp. polariton) lifetime. Wxx describes the exciton-exciton 
interaction rates. The exciton-phonon scattering rate reads: 

( E i- E k,R) 

W xp = W X pe a B 

The wave vector conservation rule is replaced by an overlap integral between the spatial 
extensions of the initial and final states. The overlap integrals change with time, together 
with the shapes and the energies of the states according to equation ([T]). The overlap 
integrals involved are defined as follows: 



hkRR = S 3 J WAk? Pr dxdy, 
IiRRR = S 3 J \i)i\ 2 p R dxdy 



(3) 



The overlap integrals obtained by summation of scattering rates with different phonon 
wavevectors read: 

hk = S I \ip*ip k \ 2 dxdy, 



These overlap integrals are dimensionless. 

Finally, the scattering rates Wxx an d Wxp describing exciton-exciton and exciton- 
phonon interaction were used as the fitting parameters of the model, together with the 
width of the gaussian spatial distribution of the reservoir which we authorize to slightly 
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FIG. 1: Schematic representation of the exciton-exciton interaction in the reservoir leading to the 
scattering towards the quantized polariton states. The reservoir has a 2D (constant) density of 
states with a thermal distribution function. 

vary above the width given by the pumping laser spot size. We used Wxx — 2 10 3 
Wxp = 10 9 w r = 3/im . 

In order to understand the terms in the Boltzmann equation, let us consider for example 
the case of the scattering from the reservoir towards a quantized state % and backward assisted 
by exciton-phonon and exciton-exciton interaction. We neglect the final state stimulation 
within the reservoir which means that we assume a weak occupation (n R (E) <C 1) of these 
states. In all cases, the energy difference between the state i and the bottom of the exciton 
reservoir is E . 

Let us first consider the most complicated case of the forward scattering assisted by 
exciton-exciton interaction (sketched in figure S-l). In order to calculate the corresponding 
scattering rate, we have to sum over all possible scattering events. The corresponding 
scattering rates reads: 

oo oo 

WxxhRRR {n l + l)^n R (E) ^ n R {E') (1 + n R (E 1 + E + Em)) « W xx I iRRR (m + 1) N 2 R 

E=0 E> 

(5) 

The rate is proportional to iVjj and is similar to a scattering rate written between individual 
states and based on particle-particle interaction. The backward process of one particle going 
from the quantized states to the reservoir reads: 
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WxxhRRRm 5^(1 + n R (E)) J2 ME' + E + E iR )(l + n R (E')) 



E=0 E'=0 



oo oo 



Wxxhunnn, ^ (1 + n R (E)) e"^^ ^ nR (E') 



E=0 E'=0 

at TT _ r iV™ TI . T SmkT 

~WxxhRRRniN R e kT kT = w xxhRRRni — j--e w = W X xUrrr „ h2 riiN R e 



n fl (0) 27r^ 2 

(6) 

The rate is similar to the one of the forward process except by the term : e fc s T /n R (0). The 
relaxation of particles toward a low energy state is therefore favored only for a minimum 
energy difference E > k R T \n(n R (0)). 

For the phonon-assisted scattering between the reservoir R and a quantized state i the 
expressions are quite similar. The scattering rates from the reservoir towards the state i 
reads: 

WxpUr (rii + nR ( E ) ( X + n P h ( E + EiR )) ~ w xpUr (th + 1) N R ( 7 ) 

E 

The scattering rates from the state i toward the reservoir reads: 

Sk T 

W XP I lR n t ^Ph (E) (1 + n R (E)) « W XP I iR n i e- E ^ hB ' 1 — |^ ( 8 ) 

_B=0 

One can note that the expressions for the scattering processes involving exciton reservoir 
and phonons are quite similar. This is explained by the fact that excitonic reservoir is 
assumed to be thermalized at the lattice temperature. 

Despite the complexity of the equations used in our simulation, we would like to underline 
that the number of fitting parameter is really limited. 

Figured presents calculations corresponding to the experiments described in fig.4 of the 
paper. The reservoir is injected in the left pillar. The difference with the calculations pre- 
sented in fig.4 is that here the scattering rates are taken constant, equal to that obtained 
at low excitation density. Thus we neglect the effect of the polariton wavefunction renor- 
malization on the scattering rates. Clearly the main features of the experiments are not 
well described, in particular condensation in the metastable state at intermediate excitation 
powers is not reproduced. 
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FIG. 2: Calculated emission pattern considering a reservoir potential proportional to the excitation 
power in the left micropillar but without taking into account the changes in the scattering rates. 

The crucial ingredient which allows catching the basic physics in our experiments is to 
take into account the overlap integrals between initial and final states in the scattering rates. 
As shown in the paper, these wave functions are strongly distorted by interactions, and this 
dramatically modifies the way polariton relaxation occurs in the system. 
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